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Abstract 

Electromagnetic wavelets are a family of 3 x 3 matrix fields W^(x') parameterized by 
complex spacetime points z = x + iy with y timelike. They are translates of a basic 
wavelet W(z) holomorphic in the future-oriented union T of the forward and backward 
tubes. Applied to a polarization vector p = — ip^, W{z) gives an anti-selfdual 
solution W{z)p derived from a selfdual Hertz potential Z(z) = —iS{z)p, where S is 
the Synge function acting as a Whittaker-like scalar Hertz potential. Resolutions of 
unity exist giving representations of sourceless electromagnetic fields as superpositions 
of wavelets. With the choice of a branch cut, S splits into a difference S'^{z) — S~{z) of 
retarded and advanced pulsed beams whose limits as y — > give the propagators of the 
wave equation. This yields a similar splitting of the wavelets and leads to their com- 
plete physical interpretation as pulsed beams absorbed and emitted by a disk source 
D{y) representing the branch cut. The choice of y determines the beam's orienta- 
tion, collimation and duration, giving beams as sharp with pulses as short as desired. 
The sources are computed as spacetime distributions of electric and magnetic dipoles 
supported on D{y). The wavelet representation of sourceless electromagnetic fields 
now splits into representations with advanced and retarded sources. These representa- 
tions are the electromagnetic counterpart of relativistic coherent-state representations 
previously derived for massive Klein-Gordon and Dirac particles. 

1 Hertz Potentials 

Acoustic and electromagnetic wavelets were introduced in [K92, K94, K94a] as generalized 
frames of localized solutions of the homogeneous wave equation and Maxwell's equations, 
giving wavelet-like representations of general sourceless solutions. They were defined in 
Fourier space, and explicit space-time expressions were found for the acoustic but not the 
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electromagnetic wavelets. Here we compute the electromagnetic wavelets in spacctime and 
show that they have a simple interpretation as pulsed beams emitted and absorbed by local- 
ized electric and magnetic dipole distributions. By separating their advanced and retarded 
parts, the wavelet representation of electromagnetic fields is generalized to include localized 
sources, interpreted as polarization distributions generating Whittaker-like Hertz potentials. 

We begin by reviewing the derivation of electromagnetic fields by Hertz potentials. Our 
presentation is essentially a translation of [N55] to the language of differential forms [AMT88, 
T96] . An electromagnetic field is represented by two 2- forms F, G satisfying Maxwell's 
equations, 

dF^O {F^E- dxdt + B-s) (1) 
dG^J (G^D-s-H-dxdt). (2) 

Here J is the current 3-form, products of forms are wedge products, dx = {dx^,dx'^,dx^), 
and s is the vector-valued spatial area form s — — * dxdt — {dx^dx^, dx^dx^, dx^dx^), where 
* is the Hodge duahty operator in Minkowski space, defined [R] so that = 

In addition to (1) and (2), we need constitutive relations, given here in terms of the electric 
and magnetic polarization densities Pg, of the medium: 

*F + G = P {P = Pe- s + Pm-dxdt), (3) 

where *F = —E-s+B-dxdt is the Hodge dual of F. Without specifying initial and boundary 
conditions, the best we can hope for is general (or local) solutions up to an arbitrary sourceless 
field. We now show to derive these using Hertz potentials. 

Since F is closed, it can be derived from a potential 1-form A subject to gauge transforma- 
tions: 

F^dA^dA', A^A^dx^", A'^A-dx- (4) 

Applying the co-differential operator 5 = *d* gives 

5 A' ^5A- Sdx ^5A- Dx, (5) 

where H — dS + 5d is the (Hodge) d'Alembertian operator on differential forms, reducing to 
the wave operator on functions. Locally at least, % can be chosen to satisfy = 6A, so 
that SA' = and A' is in the Lorenz gauge. ^ By the Poincare lemma, *A' can be derived 
from a potential 2-form: 

d*A'^0 =^ *A' ^ -dZ (Z = Ze • s + • dxdt). 

This shows that the original potential, which need not be in Lorenz gauge, can be written 
locally as 

A^*dZ + dx, F ^d* dZ. (6) 
^Due to L.V. Lorenz, not H.A. Lorentz; see Penrose and Rindler [PR84]. 
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Z is called a Hertz potential with electric and magnetic vectors Zf., Z^ ■ The relation between 
the components of F and Z can be expressed compactly in terms of the anti-selfdual form 

F~ = i * F~ and the self dual form Z+ = — i * Z+, defined by 



2F- = F + i*F = 2F-u} 2F = E + iB oj = dxdt - is 

2Z+ = Z - i * Z = 2Z • a; 2Z = Z^-iZ^ w = dxdt + is. (7) 



Namely, 



2F- = d*dZ + iSdZ = d*dZ - idSZ + iOZ 

^d*dZ-id*d*Z + iUZ ^2d* dZ+ + iUZ, 

which translates to the complex vector equation 

2F = 2iCZ - DZe (8) 

where JC is the operator 

£Z = Vx(VxZ)+ idtV X Z. 

To obtain a wave equation for Z in terms of the sources, we need a stream potential [N55] 
for J, i.e., a 2- form G* such that 

J=-dG* [G* = -D* -s + H* -dxdt). 

The resemblance ~ —G is intentional, since one possibility is G* = —G. But this would 
be circular since the field is unknown. A simple stream potential Gq can be defined as a 
time integral of the current. Then G differs from — at most by a homogeneous solution 
of (2), i.e., an arbitrary exact 2-form: 

G^-Gl + da, aeW (9) 

The constitutive relation (3) imphes 

P = *F + G = ddZ -G*o + da = DZ - d6Z - Gq + d^, 

or 



DZ = P + G* (10) 



where 

G* ^Gl-d(3, f3^a-5Z 

is another stream potential. The dependence oi P on Z is not a problem since (3, like a, is 
arbitrary. The selfdual representation of (10) is 

□Z = P + G*, 2P ^ Pm-iPe, 2G*^H* + iD*. (11) 

Thus Z is produced by a unified source composed of material polarization and stream po- 
tential; the latter may be viewed as an effective classical vacuum polarization induced by 
the current. If the medium is polarizablc, magnetizable or conducting, these sources in turn 
depend on the field. In that case it is better to use the common multiplicative form of 
the constitutive relations for the induced sources (assuming linearity), which gives a mod- 
ified version of (10). For a thorough study of Hertz potentials and their gauge theory, see 
[N55, N57]. 



3 



2 Hertz potentials in Fourier space 



Initially, the construction of electromagnetic wavelets will be based on holomorphic fields 
obtained by extending sourceless anti-selfdual solutions to complex spacetime. Eventually, 
sources will be introduced that preserve the holomorphy of the fields locally, outside the 
sources. But for now, we specialize to a homogeneous field in vacuum, 

J = p = =^ dP- = ^ idtF = V X F, V ■ F = 0. (12) 

Solutions are given in terms of selfdual Hertz potentials by (8): 

□Z(x)=0, F{x)^iCZ{x). (13) 

The Fourier solution of the wave equation is 

Z{x) ^ f dk e'^^'Zik) ^ f dk e'^''Z{k) + [ dk e^'=^Z(A;), (14) 
Jc Jc+ Jc- 

where x = (a;, i), k — {k,ko), kx — kot — k-x, C± are the positive and negative- frequency 
fight cones 

C± = {A; : ±ko = a; > 0}, u = \k\, and dk = d^k/{167r^uj) 

is the Lorentz-invariant measure on the double cone C = C+ U C_. The coefficient function 
F{k) is the restriction of the Fourier transform of F{x) to C. Inserting the definition of £, 

F{x) ^ I dk e'^'^Fik) dk 6^^=^ \-k x{kx Z)+ ik^k x z\ . (15) 
Jc Jc ^ ^ 

Note that (12) requires 

F ^inx F, where n{k) = k/ko, = 1 on C. (16) 
Thus for every k E C, F{k) is an eigenvector with eigenvalue 1 of 

E{k) : ^ defined by Ev^inxv. 

But 

E>'^v = v- n{n ■ v) ^ §^ = §, 

so § has the nondegenerate spectrum {—1, 0, 1} and the orthogonal projection to the eigen- 
space with § = 1 is 

F{k) = i + §] =^ §P = P = P* = p2. (17) 
We can now write (15) in the form 

F{x)= I dke'^'^Fik), F = 2icu¥Z = PF. (18) 
Jc 
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3 Extension to complex spacetime 



Sourceless waves are "boundary values" of holomorphic functions in a sense to be explained. 
Denote the future and past cones at a; = by 

V^^{y^{y,s)eR^:±s>\y\}. 

The forward and backward tubes [SW64, PR84] are the complex domains 

Since they are disjoint, we are free to give them independent orientations. We orient them 
oppositely, considering 7^ to have a positive and T a negative orientation. Specifically, 
let Tj_ be oriented by its volume form r = d'^x d'^y and T_ by the time-ieversed form — r 
{dx — >• dx, dt — > —dt, dy^ — >• —dy^). The oriented union (chain) 

r = r+ - r_ 

will be called the causal tube. 

We will extend sourceless fields such as F{x) to T, then interpret these extensions Fi^z) 
physically by examining their behavior on real spacetimes defined by slices 

R^y = {x + iy : X eR'^}, y e , 

in the spirit of Newman, Plebanski, Penrose, Rindler, Robinson, Trautman and others [N65, 
NJ65, N73, NW74, N02, PR78, PR84, T62]. As y 0, the extended fields converge to their 
"boundary values" on R'^ in a sense to be made specific. (It is only in this sense that F{z) is 
to be regarded an extension of F{x).) But note that R"^ is not the topological boundary of 
which would be the 7-dimenional set {x + iy : ±yo > 0, y^ — 0}, but its Shilov boundary 
[H73]. 

Solutions of homogeneous wave equations, such as Z{x) in (14), can be extended to T by 
the analytic-signal transform 

Z{z) = s [ dke'^'Z{k), z = x + iyeT, (19) 

where 

„ /C+ if5>0 
s = sgns, Cs=< {s = yo). 

I C_ if s < 0. 

Note that the positive and negative frequency parts of Z are extended to the forward and 
backward tubes, respectively. Since Z is complex to begin with, these parts are independent. 
The exponential factor e"'^^ in (19) decays as \k\ oo whether y G Kt , so if Z{k) is 
reasonable (of polynomial growth, say), then the integral defines Z{ holomorphic 
function. Since e~^y decays least rapidly along those rays k E C± that are "nearly parallel" 
to y e T4, it follows that such rays are favored by the extension (19) if y is "nearly lightlike." 
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Since Z{x) is the sum of its positive and negative frequency parts, the sign s in (19) tells us 
that it can be recovered (in a distributional sense) as the difference of the (Shilov) boundary 
values 

Z{x) ^ Z{x + ieO) - Z{x -ieO), (20) 

where 

Z{x ±iO) ^limZ{x±iey), yeV+ 

and the limit can be shown to be independent of y E V^. 
It is also possible to define Z more suggestively as 

Z{x + iy) = —j -—Z{x + Ty), (21) 

which is defined (but not holomorphic) for all y e M!^ and can be shown to reduce to (19) 
when y eV±. This definition was used in [K90, KS92, K94] witfiout the orienation s. Only 
recently have I understood the value of this orientation, which connects the transform to 
hyperfunction theory and locality (see below). 

The electromagnetic field can now be extended similarly as 

F{z) = s [ dk e'^'F{k) = iCZ{z) = 2is [ dk e*^VP(A;)Z(A;), (22) 

where C differentiates with respect to x (so as to be defined in R^). The positive and negative- 
frequency parts of F also have positive and negative tielicities [K94], so the restrictions of 
F{z) to Tj_ and T are positive and negative-hehcity solutions. 

In one dimension, the above extension is trivial: the positive and negative frequency parts 
extend to the upper and lower complex half-planes, and a quite general distribution can be 
written in the form (20). (An elementary form of this idea was applied to communication 
theory by Dennis Gabor, who called the extensions analytic signals; hence the name.) In 
dimension n > 1, we need something to replace the half-lines of positive and negative fre- 
quencies. When the support in Fourier space is contained in a double cone (like the convex 
hulls of C±), these become the appropriate replacements and the tubes over the open dual 
cones (in this case V±) replace the upper and lower half-planes. 

A hyperfunction H{x) in [K88, KS99] is a generalized distribution defined, roughly, by 
difi^crcnccs in the boundary values of a set of holomorphic functions Hk{z) in a set of complex 
domains enveloping the support of H. The functions H^^z) are called generating functions 
for H{x)? Equation (20) suggests thinking of the restrictions of Z{z) and F{z) to 7^ as 
generating functions for Z{x) and F{x). We will not attempt to make rigorous use of 
hyperfunction theory; rather, it will serve mainly as a guide. This will not only free us 

^Morc accurately, differences in boundary values must in general be replaced by sheaf cohomology. For- 
tunately, for solutions of homogeneous relativistic equations (Klein-Gordon, Dirac, etc.), simple differences 
like (20) suffice. 
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from its highly technical requirements but also allow us to go beyond it, as some of our 
constructions (like the extended delta functions S^{z) and S^{z) below) do not fit neatly into 
the theory. Also, since our goal is to understand physics directly in T, we do not view F{z) 
and Z{z) merely as tools for the analysis of the boundary distributions F{x) and Z{x), as 
they would be in hyperfunction theory. 

4 Electromagnetic wavelets 

To construct the wavelets, we need a Hilbert space of solutions. Initially the inner product 
will be defined in Fourier space. It is uniquely determined (up to a constant factor) by the 
requirement of Lorentz invariance to be 

(Fi,F2)= / ^ F{F2 = a! dku;^zlFZ2. 
Jc ^ Jc 

Note that the measure dk/cu^ oc d^k/cu^ is invariant under scaling. In fact, an equivalent 
inner product has been shown to be invariant under the conformal group C of Minkowski 
space [Gr64] . Therefore the Hilbert space of anti-selfdual solutions 

n = {F: ||Ff = {F,F) < 00} 

carries a unitary representation of C. We want to construct the wavelets so as to preserve 
the connection with the conformal group. We therefore define the wavelets in the same spirit 
as the relativistic coherent states for the free Klein-Gordon and Dirac fields [K77, K78, K87, 
K90] . For any fixed z eT, consider the evaluation map 

S^-.n^C^ defined by S^F = F{z). 

This is a bounded operator between Hilbert spaces (where is given its standard inner 
product). Its adjoint is, by definition, the electromagnetic wavelet^ 

= £* -.C^ ^n. (23) 

maps any polarization vector p G C'^ to a solution W^p e H. Therefore Wzix') must be 
a matrix-valued solution of Maxwell's equations (three columns, each a solution). But 

W*F = S^F = F{z) = s [ dk e'''^FF{k) \- FF = F, 

Jcs 

therefore, remembering the measure dk/uj"^ for the inner product in Ti, the expressions for 
in the Fourier and spacetime domains are 

W,(A;) = se-^'=VP(A;), W,{x') = s [ dk e'''^'''-'^u;''F{k). (24) 

^We arc really using a vectorial version of the Riesz representation theorem. If 7i were a Hilbert space 
of sufficiently nice scalar functions, £z would be a bounded linear functional and (1) its representation by 
the unique elelment of H guaranteed to exist by the Riesz theorem. 
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The reproducing kernel is defined as 

K{z', z) = WIM, = 9{y'y) [ dk e'''^''-'^uj'^F{k) = 9{y'y)W{z' - z), (25) 

where 9 is the Heaviside step function, the factor 9{y'y) enforces the mutual orthogonahty of 
wavelets parameterized by in the forward and backward tubes, and the holomorphic matrix 
function 



W{z) = [ dk e'^'uj^F{k) (26) 

generates the entire wavelet family by translations: 

W,{x') = sW{x' -z), zeT. (27) 
We now compute W{z) exphcitly. Applying it to a vector p e gives 

2W{z)p ^2 [ dk e'VF{k)p = ijCR{z)p, (28) 

where 



iR{z) 



I dk e'^' I dkko e'^\ ko = sgn A;o . (29) 



The second equality, which holds because s and k^ have the same sign, shows that R{z) is 
the analytic-signal transform of 



R{x) = -i I dk he'''''. (30) 
c 



It is easily checked that this integral gives the (unique) solution to the foUwing inital-value 
problem: 

DR{x) = 0, R{x, 0) = 0, dtRix, 0) = d^{x). 

R is known as the Riemann function of the wave equation [T96]. The integral is readily 
found to be the difference between the retarded and advanced propagators, 

Rix) = R+(x) - R~(x), R^ix) = (31) 

47rr 

On the other hand, we find 

S{z) = iR{z) = [ dk e''' = -T^' = ^f^^"" = - (32) 

JCs 

This is nothing but Synge's "elementary solution" [S65, p. 360] of the wave equation! Among 
other things, it has been used by Trautman [T62] to construct interesting (null, curling) 
analytic solutions of the Maxwell and linearized Einstein equations. 
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I believe this connection between the Riemann and Synge functions confirms that the trans- 
form (19) is the "right" extension of solutions to T. In the absence of the orientation s (with 
(20) now a sum), the Synge function is no longer the analytic signal of iR{x). In particular, 
we lose the connection with Huygens' principle in the limit y ^ 0. 

The Synge function is holomorphic on the complement of the complex null cone 

Af^{zeC^:z^^ 0}. 

In particular, note that 

X + iy & Af r\T =^ x^ — > 0, and xy — 0. 

But the second equality is impossible since, by the first, x and y are both timelike. Thus 
S{z) is holomorphic in T. 

According to (13), 2'W{z)p can thus be computed from the selfdual Hertz potential 

2Z{z) = R{z)p, (33) 
which gives explicit expressions to the entire wavelet family. 

• All the wavelets are translations of W{z). This can be further reduced by using the 
symmetries of W{z), which reflect those of Maxwell's equations — i.e., the conformal group. 
For example, W{z) is homogeneous of degree —4, so the Lorentz norm of y can be interpreted 
as a scale parameter: 

W,{x') = W {{x' -z)/X), A = > 0. (34) 

• It is easily shown that 

W(z)*=W(5), hence W,{x')* ^W^{x'). 

• There exist (many equivalent) resolutions of unity [K94], obtained by integrating over 
various subsets T> <zT with appropriate measures d^D'- 

[ diiTy{z)W,Wl = In. (35) 
Jv 

This is a "completeness relation" dual to the (non-) "orthogonahty" relation (25). Each 
resolution gives a representation of solutions as superpositions of wavelets, 

F{x')= [ dfij){z)W,{x')WlF ^ f dpLv{z)Wz{x')F{z), (36) 
Jv Jv 

with the analytic signal F{z) on V as the "wavelet transform." One natural subset for a 
resolution is the Euclidean region of real space and imaginary time, 

£ — {{x, is) : X E R^, s 7^ 0}, djJLsix, is) — d^x ds. 
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Since A = |s| on S, the wavelets are now parameterized by their location and scale, just as 
in standard wavelet analysis [D92], with the Euclidean time as scale parameter. 

• Applying the analytic-signal transform to (36) gives 

F{z')= [ diiviz)WlM,F{z)= [ diiT,{zMz',z)F{z), (37) 
Jv Jv 

which explains the term "reproducing kernel." 

• By construction, the wavelets transform covariantly under the Poincare group, which 
endows their parameters z with physical significance [K94]. In particular, any resolution 
(35) can be boosted, rotated or translated to give another resolution. Furthermore, since Ti 
carries a unitary representation of the conformal group C, for which T is a natural domain, 
it is useful to study the action of C on wavelets. To some extent this has been done in [K94]. 
A connection has been discovered recently with work on twistor-like transforms by Iwo 
Bilynicki-Birula [B02], and we have begun a joint project together with Simonetta Prittelli. 

5 Interpretation as Hertzian pulsed beams 

Recall that the absence of sources was the price we paid at the outset for the holomorphy used 
to construct wavelets. But now that we have the wavelets, it turns out that we can make 
them physical (causal!) by introducing sources in a way that preserves their holomorphy 
everywhere outside these sources - which is the most that can be expected! 

The introduction of sources will be based on a causal splitting of R{z) similar to that of 
R{x). But whereas the splitting of R{x) is Lorentz-invariant [the hyperplane {t = 0} can be 
tilted without affecting R^{x)], the corresponding splitting of R{z) is frame-dependent. We 
must choose a 4- velocity v {v^ = l,Vo > 0) and split z into orthogonal temporal and spatial 
parts 

I 2 2 2 2,2 

Z = TV + Zs, T = VZ, Z = Zq — Z = T -\- Zg . 

Without loss of generahty, we take v — (0, 1) in the following. The general case is recovered 
by letting 

2 2 2 

2^0 —> T, Z — > T — Z . 

Begin with the factorization 

z^-z^ = (zo-f)(zo + f), 

where 

f{z) — -s/z^ — r"^ — a? -\- 2ix ■ y, z — x + iy, r — \x\, a — \y\. 

For motivation, think of f as the complex distance from a "point source" at —iy to the 
observer at x. A study of potential theory based on complex distances in C", and the 
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connections they furnish between elhptic and hyperboUc equations, is in progress; see [KOO]; 
see also the related papers [KOI, KOla, N02]. Given y ^ 0, the branch points of f form a 
circle in the plane orthogonal to y, 

C{y) = {a; e : r = a, a; • y = 0}, 

and following a loop that threads C changes the sign of f. We make f single-valued by 
choosing the "physical" branch defined by 

Re f > 0, so that t/ ^ =^ f ^ +r. 




Complex distance: The real part (left) and imaginary part (right) of f{x+iy) as functions 
oi X = (xi, 0,2:3), with y = (0,0, 1). The real part is a pinched cone, with the branch disk 
in the X1-X2 plane projected to the interval [—1, 1] of the xi axis. 

The resulting branch cut is a disk^ bounded by C, 

D{y) = {x : r < a, x ■ y = 0}. 

The extended spatial delta function is defined as the distributional Laplacian with respect 
to X of the holomorphic Coulomb potential: 

6%z) = -A-^. (38) 

As a distribution in x, it is a natural extension of the point source S^{x). Roughly speaking, 
displacing the singularity to —iy opens up a "light cone" in with the y-axis as "time," 
of which the branch circle is a "wave front" and D its interior. Similar remarks apply to 
extended delta functions in C", which have provided an intriguing connection [KOO] between 
fundamental solutions of Laplace's equation in Eulidean M" and the initial- value problem 
for wave equations in Lorentzian R". 

^An equivalent branch cut is obtained by continuously deforming the disk to an arbitrary membrane M 
bounded by C. 
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Returning to spacetime, R{z) decomposes into partial fractions: 



^i^) = T^T ^ ^ ^^(^) - ^"(^)' (39) 



where 



Ri(z)^^-1 -. (40) 

Note that although R{z) = —iS{z) is holomorphic in T, its retarded and advanced parts 
R^{z) are not. Viewed on the slice , they are singular on the world-tube D traced out by 
the branch cut D. Everywhere in T outside this source region, they are still holomorphic. 

The boundary values of are given by the Plemelj formulas 

R^{x + zO) = — = ' + V- 



where V denotes the Cauchy principal value. Therefore the jumps across are 

R^(x + iO)-R^(x-iO) ^ ^^[^''^ ^R^(x). (41) 

47rr 



We now show that R'^iz) have very interesting physical interpretations even when y 7^ 0, by 
looking at their behavior in slices . 

Guided by the successful definition (38) of the extended spatial source, define the extended 
spacetime delta function 

S\z) = DR^{z), (42) 

where □ is the distributional d'Alembertian acting on x. Since — OR' = DR = 0, 
there is no sign ambiguity on the left. A detailed study of S'^{x + iy) is somewhat involved 
and will be given elsewhere [K03] . Here we note the following. 

• 5^(x + iy) is a well-defined Schwartz distribution in R^. 

• By (41), 

5^ {x + iO) - (x - iO) = Di?* (x) = 5^ (x) . (43) 



• It is easy to show that □-R*(x + iy) = at all points of regularity. Thus 6^{x + iy) is 
supported on the world tube D representing the evolution of the source disk in R^. 

• In spite of (43), the restrictions of 5^ (2;) to 7^ are not generating functions for (5'^(a;) since 
they are not holomorphic in any neighborhood of R"^. (This is what 1 meant earlier by saying 
that not all our constructions fit neatly into hyperf unction theory.) 
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• In the far zone, we have 

r ^ a f — y/r"^ — a? -\- 2iar cos 6 ~ r + ia cos 6 
{zQ — t — is, \x\ — r, \y\ — a, x-y — racos9), 

giving a simple expression for the far field from which the pulsed-beam interpretation can 
be read off easily: 



r ^ a 



STT^r t^r -\-i{s^ a cos 9) 






Time-lapse plots at i = 1, 5, 10, 15, 20 of \R^{x + iy)\ with y — (0, 0, a, 1) as functions of 
X — {xi, 0, X3). Clockwise from upper left: a — .5, .9, .99, .999. The pulses propagate along 
y — (0, 0, a), getting sharper and sharper as a — > s = 1. 

• R'{x + iy) is an advanced pulsed beam converging toward D along y/s and absorbed in 
D around t = 0. 

• R^{x + iy) is a retarded pulsed beam emitted from D around t = and propagating along 

• D acts like an antenna dish, simultaneously absorbing and emitting resulting in 
the sourceless pulsed beam R{x + iy) focused ai x — 0. 
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• Both beams have duration \s\—a along the beam axis and are focused as sharply as desired 
by letting (y/s)^ 1. As {y/s)"^ — > 0, they become spherical pulses of duration 

Inserting (40) into (33) and (28) gives 

W(z) = W+(^) - W-(z), 

where 

{z)p = iCZ^ (z) , 2Z* (z) = Ft {z)p. (44) 
By (42), the polarization associated with this Hertz potential is 

2P{z) = 2DZ^{z) = p5\z). (45) 
By (43), the polarization on is an impulsive dipole 

Pm{x) - iPe{x) = 2P{x) = 2P{x + iO) - 2P{x - iO) =p6'^{x), 
giving an intepretation of p in terms of magnetic and electric dipole moments 

P^Pm- iPe- (46) 

Of course, the extension P{x) P{z) mixes the electric and magnetic dipoles since S'^{z) is 
complex. Nevertheless, we are free to define real polarizations and causal/ anticausal Hertz 
potentials and fields in the slice by 

PUz)-tPe{z)=p~S\z) 

ZUz)-iZt{z)=R^{z)p 
E^{z) + iB^{z) = iC[R^{z)p]. 

The left-hand sides are then solutions of Maxwell's equations inheriting the pulsed-beam 
interpretations derived for R'^{z), but now in the context of electrodynamics instead of 
scalar wave equations. 

Whittaker [W04, N55] has shown that given any constant, nonzero vector v e R^, a general 
Hertz potential can be reduced to two scalar potentials He , H^ in the form 

Ze{x) = Ile{x)v and Zm{x) = Il.rn{x)v, Or 2Z{x) = {lirn — iTl.e)v. 

Comparison with (44) shows that R^{z) are scalar Hertz potentials of Whittaker type for 
2W±(z)p. 

Finally, the wavelet decomposition (36) of sourceless anti-selfdual fields splits into 

F{x') = F+{x') - F-(x'), F^(x') = / diiv(z)Wt(x')F(z). 

Jv 

F^ and F' are fields absorbed and emitted, respectively, by source disks distributed in T) 
according to the coefficient function F[z). This may be used for analyzing and synthesizing 
fields with general sources. 

Pulsed beams similar to the above have also been studied and applied in the engineering 
literature, from a different point of view; see [HFOl] for a recent review. 
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